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Monomials function has always been considered as a significant and most extensively used function in 
real living. Resource allocation, structure optimization and technology management can often apply 
these functions. In optimization problems the objective functions can be considered by monomials. In 
this paper, we present monomials geometric programming with fuzzy relation inequalities constraint 
with max-product composition. Simplification operations have been given to accelerate the resolution 
of the problem by removing the components having no effect on the solution process. Also, an algorithm 
and a few practical examples are presented to abbreviate and illustrate the steps of the problem 
resolution. 

© 2008 Elsevier Ltd. All rights reserved. 



1. Introduction 

Fuzzy relation equations (FRE), fuzzy relation inequalities (FRI), 
and their connected problems have been investigated by many 
researchers in both theoretical and applied areas (Czogala, Drew- 
niak, & Pedrycz, 1982; Czogala & Predrycz, 1982; Di Nola, Pedrycz, 
& Sessa, 1984; Fang & Puthenpura, 1993; Guo, Wang, Di Nola, & 
Sessa, 1988; Gupta & Qi, 1991; Higashi & Klir, 1984; Han, Song, & 
Sekiguchi, 1995; Hu, 1998; Li & Fang, 1996; Perfilieva & Novak, 
2007; Prevot, 1985; Sheue Shieh, 2007; Wang, 1984, 1991; Zadeh, 
1965; Zener, 1971). Sanchez (1977) started the development of the 
theory and applications of FRE treated as a formalized model for 
non-precise concepts. Generally, FRE and FRI have a number of 
properties that make them suitable for formulizing the uncertain 
information upon which many applied concepts are usually based. 
The application of FRE and FRI can be seen in many areas, for in- 
stance, fuzzy control, fuzzy decision making, systems analysis, fuz- 
zy modeling, fuzzy arithmetic, fuzzy symptom diagnosis and 
especially fuzzy medical diagnosis, and so on (see Adlassnig, 
1986; Berrached et al., 2002; Czogala and Predrycz, 1982; Di Nola 
and Russo, 2007; Dubois and Prade, 1980; Loia and Sessa, 2005; 
Nobuhara et al., 2006; Pedrycz, 1981, 1985; Perfilieva and Novak, 
2007; Vasantha Kandasamy & Smarandache, 2004, chap. 2; Wen- 
stop, 1976; Zener, 1971). 

An interesting extensively investigated kind of such problems 
is the optimization of the objective functions on the region 
whose set of feasible solutions have been defined as FRE or FRI 
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constraints (Brouke & Fisher, 1998; Fang & Li, 1999; Fernandez 
& Gil, 2004; Guo & Xia, 2006; Guu & Wu, 2002; Higashi & Klir, 
1984; Khorram & Hassanzadeh, 2008; Loetamonphong & Fang, 
2001; Loetamonphong et al., 2002; Zadeh, 2005). Fang and Li 
solved the linear optimization problem with respect to FRE con- 
straints by considering the max-min composition (Fang & Li, 
1999). The max-min composition is commonly used when a sys- 
tem requires conservative solutions, in the sense that the good- 
ness of one value cannot compensate for the badness of another 
value (Loetamonphong & Fang, 2001). Recent results in the liter- 
ature, however, show that the min operator is not always the 
best choice for the intersection operation. Instead, the max- 
product composition has provided results better than or equiva- 
lent to the max-min composition in some applications (Adlass- 
nig, 1986). 

The fundamental result for fuzzy relation equations with max- 
product composition goes back to Pedrycz (1985). A recent study 
in this regard can be found in Brouke and Fisher (1998). They ex- 
tended the study of an inverse solution of a system of fuzzy rela- 
tion equations with max-product composition. They provided 
theoretical results for determining the complete sets of solutions 
as well as the conditions for the existence of resolutions. Their re- 
sults showed that such complete sets of solutions can be character- 
ized by one maximum solution and a number of minimal solutions. 
An optimization problem was studied by Loetamonfong and Fang 
with max-product composition (Loetamonphong & Fang, 2001) 
which was improved by Guu and Wu by shrinking the search re- 
gion (Guu & Wu, 2002). The linear objective optimization problem 
with FRI was investigated by Zhang, Dong, and Ren (2003), where 
the fuzzy operator is considered as the max-min composition. 
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Also, Guo and Xia presented an algorithm to accelerate the resolu- 
tion of this problem (Guo & Xia, 2006). 

The geometric programming (GP) theory proposed in 1961 by 
Zeneret al. for first time (Duffin, Peterson, & Zener, 1967; Peter- 
son, 1976). Business administration, economic analysis, resource 
allocation and environmental engineering have a large number 
of applications in GP (Zener, 1971). The fuzzy geometric pro- 
gramming problem proposed by Cao (2001). He considered a 
few number of power systems problems (Cao, 1999) and also 
Liu applied it in economic management (Liu, 2004). The fuzzy 
geometric programming with multi-objective functions has stud- 
ied by Biswal (1992) and Verma (1990). In order to show impor- 
tance of geometric programming and the fuzzy relation equation 
in theory and applications a fuzzy relation geometric program- 
ming problem has proposed by Yang and Cao (2005a, 2005b). 
Furthermore, they discussed optimal solutions with two kinds 
of objective functions based on the fuzzy max-product operator. 
Also, they consider monomial geometric programming with fuz- 
zy relation equation (FRI) constraints with max-min composition 
(Yang & Cao, 2007). 

In this paper, we consider the monomial geometric program- 
ming of the FRI with the max-product operator. This problem can 
be formulated as follows: 

n 

min c 

j = i 

s.t. A x ^ d 1 (1) 

Bx^d 2 

X G [0, l] n 

where c,o, g R,c > 0,A = (a») mx „,a» e [0,1], B = (by) lxn ,b 9 e [0,1], 
are fuzzy matrices, d 1 = (d?) mxl g [0, l] m , d 2 = (d 2 ) lxl g [0, l] 1 are fuz- 
zy vectors, x = (X/) nxl g [0, l] n 0 is an unknown vector, and “’’denotes 
the fuzzy max-product operator as defined below. Problem (1 ) can be 
rewritten as the following problem in detail: 

n 

min cf]*? 

j = i 

s.t. di -x ^ d) i'gZ 1 = {1,2,..., m} ( 2 ) 

bi-x^d 2 iel 2 = { 1,2,. ..,4 
O^Xj^l j ej = {1,2,..., n} 

where a, and b t are the ith row of the matrices A and B, respectively, 
and the constraints are expressed by the max-product operator def- 
inition as: 

di-x = max{ay • xf) ^ d) Vi g Z 1 

jeJ ( 3 ) 

bi • x = max{hy • xf} ^ df Vi el 2 

In Section 2, the set of the feasible solutions of Problem (2) 
and its properties are studied. A necessary condition and a suffi- 
cient condition are given to realize the feasibility of Problem (2). 
In Section 3, some simplification operations are presented to 
accelerate the resolution process. Then, in Section 4 an algorithm 
is introduced to solve the problem by using the results of the pre- 
vious sections, and a numerical example is also given to illustrate 
the algorithm in this section. Finally, the conclusion is stated in 
Section 5. 

2. The characteristics of the set of feasible solution 

Notation. We shall use, during the paper, these notations as 
follows: 



S(A,d 1 ) l - = |xg [0, l] n : cii x ^ d- j for each i gZ 1 
S(B , d 2 ) i = jx g [0, l] n : bi • x ^ d 2 jfor each i g Z 2 
S(A,d 1 ) = n S(A,d 1 ) I - = {x e [0,l] n :A-x ^ d 1 } 

iel 1 

S(B, d 2 ) = n S(B, d 2 \ = {x g [0, l] n : B ■ x < d 2 } 

ieT 

S(A,B,d 1 ,d 2 )=S(A,d 1 )nS(B,d 2 ) = {xG[0,l] n :Ax^d 1 ,Bx<d 2 }. 

Corollary 1. x g S(A,d 1 ) i for edch i e Z 1 if and only if there exists 
some j t g J such thdt Xj. ^ similorly, x e S(B , d 2 ), for edch i e 1 2 if 
and only if Xj ^ Vj e J. 

Proof. This clearly results from relations (3). □ 

Lemma 1. 

(a) S(A,d 1 )^</) ifond only if for edch i e Z 1 there exists somej t e J 
such thdt dtj. ^ d). 

(b) IfS(A,d') then T = [1,1,..., l] t lxn is the gredtest element in 
set S(A, d 1 ). 

Proof. 

(a) Suppose S(A , d 1 )#0 and x g S(A, d 1 ). Thus, x g S(A , d 1 ) f , Vi g Z 1 
and then for each i e Z 1 we have Xj. ^ for some j t g J from 
Corollary (1). Therefore, since x g S(A,d 1 ) thenx g [0, l] n and 
hence < 1, Vi g Z 1 , which implies that there is a j, g J such 

that d[j { ^ d ) , Vi g Z 1 . Conversely, suppose that there exists 
some j, g J such that dy. ^ d ] , Vi g Z 1 . Set 

x = T = [1,1,..., l]i xn , since XG[0,l] n and x j . = l>‘t, 

Vi g Z 1 then x g S(A, d 1 ),, Vi g Z 1 , from Corollary (1), and thus 
x g S(A,d 1 ). 

(b) Proof is attained from part (a) and Corollary (1). □ 

Lemma 2. 

(a) S(B, d 2 )^</>. 

(b) The smdllest element in set S(B , d 2 ) is 0 = [0, 0, ... , 0] t lxn . 

Proof. Set x = 0 = [0, 0, ... , 0] t lxn . Since d 2 ^ 0 and by ^ 0 (in case 
bi = 0 2 the problem is always well-defined and it is clear), 
then|- ^ 0. Therefore x, < ^r,Vi g Z 2 ,j g J, then Corollary (1) 
implies that x g S(B, d 2 ) and hence parts (a) and (b) are proved. □ 

Theorem 1 ( Necessnry condition ). If S(A,B,d\d 2 )^(/), then for edch 
i g Z 1 there exist j e J such that dy ^ d] 

Proof. Suppose that S(A,B, d 1 , d 2 )^</>, then since 
S(A,B, d\d 2 ) = S(A,d 1 ) nS(B, d 2 ), therefore S(A,dV</>. Now, the 
theorem is proved by using part (a) of Lemma (1 ). □ 

Definition 1. Set x = (x,) nxl where 

1 Vi : by ^ d 2 

min : by > d 2 j otherwise 

Lemma 3. Z/S(B, d 2 )^</> tZien x is the gredtest element in set S(B , d 2 ). 
Proof. See [19, p. 348]. □ 
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Corollary 2. S(B, d 2 ) = {x e [0, l] n : B ■ x ^ d 2 } = [0,x], in which x Lemma 5. Suppose S(A,d 1 )^(/) then S(A,d 1 ) = u[x(e),l] where 
and 0 are as defined in Definition ( 1 ) and Lemma (2), respectively. X(e) = (x(e) : e e j/}. A e 



Proof. Since S(B, d 2 )^</>, then 0 and x are the single smallest ele- 
ment and greatest element, respectively, from Lemmas (2) and 
(3). Let x g [0,x], then x g [0, l] n and x ^ x, Thus, b { ■ x < b, • x ^ d 2 , 
Vi g I 2 which implies x e S(B, d 2 ). Conversely, let x e S(B , d 2 ) from 
part (b) of Lemma (2), 0 ^ x and also x g S(B , d 2 ) f , Vi g I 2 . Then, Cor- 
ollary (1) requires Xj ^ Vi g / 2 and Vj g J. Hence, Xj ^ x,, Vj g J , 
which means x ^ x. Therefore x g [0,x]. □ 



Definition 2. Let j, = {j g J : a*,- ^ d-}, Vi g I 1 . For each j e j„ we 
define i x(/) = (i x(/)fc ) nxl such that 



**(/)* ” I ^ k j 

0 k#j 



Lemma 4. Consider a fixed i el 1 . 

(a) If d 1 then the vectors i x(/) are the only minimal elements of 
S{A : d l ) i for eachjefi. 

(b) If d 1 = 0 then 0 is the smallest element in S(A , d 1 ),-. 



Proof. 

(a) Suppose j e j, and i g J 1 . Since i x(j) . = then i x(/) g S^d 1 /, 
from Corollary (1). By contradiction, suppose XGS^d 1 ^ 
and x < i X 0). Hence we must have x, < ^ and x k = 0 for 

kej and fc^j. Then x, < Vj g j, and hence x $ S(A , d 1 ), from 

Corollary (1), which is a contradiction. 

(b) It is clear from Corollary ( 1 ) and the fact thatx, ^ 0, Vj g j. □ 

Corollary 3. If then S(A,d 1 ) i = {xg [0, l] n : a, x ^ 

d/} = Uj e j.[i x (j), 1], where i g J 1 and i x(/) is as defined in Definition (2). 

Proof. If S(A,d 1 ) i ^(j) then from Lemmas (1) and (4), vector 1 is the 
maximum solution and the vectors i x( j } , Vj g J t are the minimal 
solutions in S(A,d l ) i . Let xg u [i x( j). 1). Then xg [i x ^,l] for some 
j g fi and thus x g [0, l] n and% ^ i x(j) . = |; from Definition (2), 
hence, x g S(A,d 1 ) i from Corollary (1). Conversely, let x g S(A,d 1 ) i . 
Then there exits some / g j such that Xy ^ ^ from Corollary (1). 
Since xg [0,lf, then £-<1, and thus / g j,. Therefore, 
i X (/') ^ x ^ 1 which implies x g u [i X (j), !]• □ 

Definition 3. Let e = (e(l), e(2), . . . e(m)) ej 1 xf 2 x---xj m such 
that e(i) = j g j,. We define x(e) = (x(e)j) nxl , in which x(e/ = 

max{i X (e(f)).} = max if and x(e/ ; = 0 if IJ = 0, where 

Ij = {i e I 1 ; e(i) =jj. 



Corollary 4. 

(a) If d) = 0 /or some i el 1 , then we can remove the ith row of 
matrix A with no effect on the calculation of the vectors x(e) 
for each e ef =} x x j 2 x ■ ■ ■ x J m . 

(b) Ifj f J if Vi g /\ then we can remove the jth column of matrix A 
before calculating the vectors x(e), Ve g ] It and set x(e) j = 0/or 
each e g j/. 

Proof. 



(a) It is proved from Definition (3) and part (b) of Lemma (4), 
because we will get the minimal elements of S(A,d 1 ). 

(b) It is proved by only using Definition (3). □ 



Proof. If S(A,d 1 )^4>, then S(A,d 1 ) i ^4>, Vi g I 1 . Therefore, we have 

S(A,d 1 )= n S(A,d 1 ) i = n [u[i X 0),T]] = n [ u [(*«», T]] 

iel 1 iel 1 jeii iel 1 e(i)eJi 

= u [ n [i x( e(0), 1]] = u [max{i x(e(i)) }, 1] = u [x(e), T] 

ecJi iel 1 e £Ji iel 1 e ^Ji 

= U [x(e),l] 

X(e) L 

from Corollary (3) and Definition (3). 

From Lemma (5), it is obvious that S(A,d 1 ) = u [x(e),l] and 

i i X 0 (e) 

X 0 (e) = S 0 (A,d ), where X 0 (e) and S 0 (A,d ) are the set of minimal 
solutions in X(e) and S(A, d 1 ), respectively. □ 

Theorem 2. //S(A,B, d 1 , d 2 )^(/), then S(A,B,d\d 2 ) = u [x(e),x]. 

X 0 (e) 

Proof. By using Corollary (2) and the result of (5), we have 

S(A,B,d\d 2 ) = S(A,d') nS(B,d 2 ) = { u [x(e), T]1 n f0,x] 

L x o(e) J 

= u [x(e),x] 

X 0 (e) 

and the proof is complete. □ 

Corollary 5 ( Necessary and sufficient condition ). S (A, B , d 1 , d 2 )^4> if 
and only if x e S^A^d 1 ). Equivalently, S(A,B, d 1 , d 2 )^(j) if and only if 
there exists some e g j 7 such that x(e) ^ x. 

Proof. Suppose that S(A,B, d/d 2 )^^, then S(A,B, d\d 2 ) = 
u [x(e),x] by Theorem (2), thus x e S(A,B,d 1 ,d 2 ), and hence 

X 0 (e) 

xeS(A,d 1 ). Conversely let xeS(A,d 1 ). Meanwhile we know 
x € S(B. d 2 ), therefore x e S(A, d 1 ) n S(B, d 2 ) = S(A, B,d\d 2 ). □ 

3. Simplification operations and the resolution algorithm 

In order to solve Problem (1), we first convert it into the two 



sub-problems below: 




min c n *j j 

jeR + 

s.t. A x = b 


(4a) 


XG [0,l] n 




min c n x ? 

i&R 

s.t. A x = b 


(4b) 



XG [0,l] n 

where R + = {j\ 0 Cj ^ 0 ,j g j} and R~ = {focj < 0 ,j g j}. 



Lemma 6. The optimal solution of Problem (4b ) is x. 

Proof. In objective function (4b) oq < 0 therefore, xj j is a monotone 
decreasing function of Xj in the interval 0 < Xj < 1 for each j g R~. 
As a result jeR -Xj j is too. Hence, x is the optimal solution because 
x is the greatest element in set S(A, B, d 1 , d 2 ). □ 

Lemma 7. The optimal solution of Problem (4a) belongs to X 0 (e). 

Proof. In objective function (4a), oq ^ 0; therefore, xj j is a mono- 
tone increasing function of x, in the interval 0 < Xj < 1 for each 
j g R + . As a result Yl jeR + Xj j is too. Now, suppose that y g S(A, 
B,d\d 2 ) is selected arbitrarily then there exists x(e 0 ) GX 0 (e) such 
that y ^ x(e 0 ). Since Yl jeR +Xj j is a monotone increasing function 
of Xj, then Yl jeR +yJ ] > Y[j e R+ x { e o)J ] \ therefore, one of the elements 
of X 0 (e) is the optimal solution of Problem (4a). □ 
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Theorem 3. Assume that x(e 0 ) is an optimal solution (not necessary 
unique) of Problem (4a), then the optimal solution of Problem (\) is 
x * defined as follows : 

j e R ~ 

1 \ x(£o)j J e R 

Proof. Consider S(A, B, d 1 , d 2 ), then by Lemmas (6) and (7) we have 

nv= n ’?■>?> n 

jej keR + ,leR~ keR + ,leR~ jej 

Therefore, x* is the optimal solution of Problem (1 ) and the proof is 
completed. 

For calculating x* it is sufficient to find x and x(e 0 ) from 
Theorem (3). While x is easily attained by Definition (1), x(e 0 ) is 
usually hard to find. Since X 0 (e) is attained by pairwise comparison 
between the members of set X(e), then the finding process of set 
X 0 (e) is time-consuming if X(e) has many members. Therefore, a 
simplification operation can accelerate the resolution of Problem 
(4a) by removing the vectors e e / such that x(e) is not optimal in 
(4a). One of such operations is given by Corollary (4). Other 
operations are attained by the theorems below. □ 

Theorem 4. The set of feasible solutions for Problem (\), namely 
S(A,B,d\d 2 / is non-empty if and only if for each i e Z 1 set 
Ji = { j */} is non-empty, where x is defined by Definition 

(n 

Proof. Suppose S(A,B,d\d 2 )^4>. From Corollary (5), 
x e S(A,B,d\d 2 ) and then we have x e S(A,d 1 ) i1 Vi e Z 1 . Thus, for 
each i e Z 1 there exists some j ej such that x, ^ from Corollary 

J Ufj 

(1), which means Vi e Z 1 . Conversely, suppose/^, Vi e l\ 
Then there exists some j ej such that Xj ^ §;, Vi e Z 1 . Hence, 
x e S(A,d 1 ) i , Vi e Z 1 from Corollary (1), which implies x e S(A,d l ). 
These facts together with Lemma (3) imply x e S(A,B,d\d 2 ), and 
therefore S(A,B, d\ d 2 )^4>. □ 

Theorem 5. If S(A,B,d* ,d 2 )^</>, then 

S(A,B, d 1 ,d 2 ) = u [x(e),x] where 

x(e) = {x(e) : e eji*=ji x J 2 x • • • x j m }. 



Proof. By considering Theorem (2), it is sufficient to show 
x(e) $ S(A,B,d\d 2 ) if e<£j/. Suppose e$ji. Thus, there exist i! el 1 
and / e J? such that e(i') =f and Xy < Then f e Ty and by Defi- 
nition (3) we have x(e).* = maxT^-1 ^ p- > x,/. Therefore, 
_ J i^f l a if J a i'f Y 2 

x(e) ^ x is not correct, which implied x(e) $ S(A, B,d ,d ) by Theo- 
rem (2). 

From the defined notation of Theorem (4), /• cf, Vi e Z 1 , which 
requires X(e) cX(e). Also, S 0 (A,B,d\d 2 ) cX(e) by Theorem (4), in 
which S 0 (A,B,d\d 2 ) is the set of the minimal elements of 
S(A,B,d\d 2 ), thus Theorem (5) reduces the search region to find 
set So(A,B,d\d 2 ). □ 

Definition 4. Let j\,j 2 e J , ay > 0 and a j2 > 0. j 2 is said to dominate 
/ if and only if 

(a) j, e J, implies j 2 e j„ Vi e l\ 

(b) For each i e Z 1 we have ^ (^-j , such that/ e Ji. 



Theorem 6. Suppose thatj 2 dominates / forj 1 ,j 2 e J, then, the min- 
imum value of the objective function is zero. 



Proof. Suppose x(eo) is the optimal solution in (4a), define 
e' = (e'(i)) mxl such that 



e'(i) = 



e 0 (0 i $ //“ 



h 



i e I-°. 
j i 



It is obvious that I? = </> and then x(e') j = 0; also, it is feasible. 
Since Yl jeJ xJ ] ^ 0 for" each x e S(A,B,d\d 2 ), and Y\ jeJ x(e')J j = 0, 



thereforex(e') is an optimal solution and the minimum value of 
the objective function is zero. □ 



4. Algorithm for finding an optimal solution and examples 



Definition 5. Consider Problem (1). We call A = (aij) mxn and 
B = (B«)i v „ the characteristic matrices of matrix A and matrix B, 
respectively, where ay = ^ for each i e I and j e J, also by = ^ for 
each i e I 2 and j e J. (Set 1 and § = oc). 



Algorithm. Given Problem (2), 

(1) Find matrices A and B by Definition (5). 

(2) If there exists i e Z 1 such that dy > 1, Vj ej, then stop. Prob- 
lem (2) is infeasible (see Theorem (1)). 

(3) Calculate x from B by Definition (1). 

(4) If there exists i e Z 1 such that d] = 0, then remove the i’th 
row of matrix A (see part (a) of Corollary (4)). 

(5) If a,) > Xj, then set ay = 0, Vi e Z 1 and Vj e J. 

(6) If there exists i e Z 1 such that dy = 0, Vj e J , then stop. Prob- 
lem (2) is infeasible (see Theorems (4) and (5)) 

(7) If there exists / e J such thatdy/ = 0, Vi e Z 1 , then remove the 
jth column of matrix A (see part (b) of Corollary (4)) and set 
x(e 0 )y = 0. If/ e R + then Ve e J If x(e) is the optimal solution 
of (4a) and the minimum value of the objective function is 
zero. Then stop. 

(8) If j 2 dominates /, (/ ,j 2 e R + ) then remove column / 
from A, Vji,j 2 ej (see Theorem (6)), and set x(e 0 ) J - 1 =0 
then Ve e J h x(e) is the optimal solution of (4a) and 
the minimum value of the objective function is zero. 
Then, stop. 

(9) Let jr = {j e j : dyV0}and/ new =JT xj 2 new x • • • xj n ™_. Find 
the vectors x(e), Ve e / new , by Definition (3) from A, and 
x(e 0 ) by pairwise comparison between the vectors 

x(e). 

(10) Find x* from Theorem (3). 



Example 1. Consider the problem below: 



minZ 


= (*i)~ 


' 2 (*2)(*3 


0 2 (*4) 


1 

2 








- 0.5 


0.8 


0.35 


0.25" 




"xr 




- 0.4 - 


0.9 


0.92 


0.9 


1 








0.9 


0.2 


1 


0.45 


0.4 




^3 




0.8 


.0.55 


0.6 


0.8 


0.64. 




_X 4 _ 




.0.65. 



"0.6 0.5 0.1 O.r 
0.2 0.6 0.6 0.5 
_0.5 0.9 0.8 0.4_ 



"xr 




-0.48- 

0.56 






_0.72_ 


_x 4 _ 







O^Xj^l, j = 1,2, 3, 4 



Step I: Matrices A, B are as follows: 
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0.8 


0.5 


1.14 


1.6 " 




1 


0.97 


1 


0.9 




4 


0.8 


1.77 


2 




.1.18 


1.08 


0.81 


1.01 _ 




‘ 0.8 


0.96 


4.8 


4.8 ' 




2.8 


0.93 


0.93 


1.12 




1.44 


0.8 


0.9 


1.8 




2-3: 


Vector x is as 


follows: 



x = [0.8, 0.8, 0.9,1] 

Step 4-5: By this step, matrix A is converted to the following: 
"0.8 0.5 0 0 ' 

- 0 0 0 0.9 

0 0.8 0 0 

_ 0 0 0.81 0 . 

Step 6-9: J n ™ = {1 , 2 }J n 2 ew = {4 }J n ™ = {2} and 

jnew = {3} Hence> x 0 (e) = {(0 0.8 0.81 0.9)} and 

x(e 0 ), the optimal solution of Problem (4a), is clearly 
(0 0.8 0.81 0.9). 

Step 10: The optimal solution of the problem is 
x * = (0.8 0.8 0.81 1) and the minimum value of 

the objective function is Z = (0.8)“ 2 (0.8)(0.81)2 
(1)“*= 1.125. 



Example 2. To show high efficiency of the presented algorithm 
in this example, we consider the problem with fuzzy relation 
equations (FRI) constraints as follows (see Example 5.1 in 
Yang & Cao, 2007 that has been solved with max-min 
composition): 



min Z = X1X2X3X4X5 
s.t. 





"0.3391 


0.4757 


0.4403 


0.5857 


0.4329" 




"xr 




0.3682 


0.3823 


0.6001 


0.4295 


0.1488 




*2 


x = 














*3 




0.6702 


0.9954 


0.6981 


0.4027 


0.8493 








0.8195 


0.6934 


0.3742 


0.0096 


0.7798 




X 4 












-*5- 



"0.5456" 

0.5244 _ 

0.8987 “ 

0.7544 

The above constraint Ax = b is equivalent with A x and 
A • x ^ b, therefore for this problem A = B, and then we can use 
the above algorithm. 



Step 5: 

- 0 0 0 0.93 0 - 

0 0 0.87 0 0 

0 0.90 0 0 0 

.0.92 0 0 0 0.96. 

Step 6: 

Step 7: 

Step 8: 

In according to this step first column dominates fifth column, 
therefore x(e 0 ) 5 = 0, and optimal solution is x = (0.92 0.90 
0.870.930) and the minimum value of the objective function is 
zero. 

We note that Lu and Fang (2001) obtained the objective value 
this example / = 0.00000011225 using a genetic algorithm and 
after 2004 iterations ( Yang & Cao, 2007). Further; Example 2 
shows the presented algorithm achieve the solution easier with 
respect to that of Yang and Cao (2007 ).Therefore, the existing ap- 
proach is superior with respect to other earlier works which has 
done so far. 



Example 3. Consider another example as follows: 



min Z — Xl *2*4*6 

s.t. 



- 1 


0.2 


0.5 


0.5 


0.4 


0.2- 




"*r 




-O.4- 


0.4 


0.2 


0.2 


0.4 


0.5 


0.5 




*2 




0.5 


0.1 


0.8 


0.3 


0.6 


0.3 


0.2 




*3 




0.3 


0.3 


0.4 


0.5 


0.8 


0.8 


0.6 




x 4 




0.6 


0.3 


0.2 


0.4 


0.5 


0.4 


0.3 




*5 




0.3 


_0.5 


0.6 


0.3 


0.3 


0.5 


0.4_ 




_*6_ 




.0.6. 


-0.3 


0.2 


0.4 


0.5 


0.5 


0.3- 






-O.3- 


0.5 


0.6 


0.3 


0.3 


0.5 


0.4 




"xr 




0.6 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 




*2 




0.6 


0.6 


0.5 


0.4 


0.3 


0.2 


0.1 




*3 


< 


0.6 


0.4 


0.4 


0.4 


0.3 


0.3 


0.3 




X 4 




0.4 


0.5 


0.5 


0.6 


0.6 


0.7 


0.7 




*5 




0.7 


0.8 


0.8 


1 


1 


1 


0.8 




_*6_ 




0.8 


.0.3 


0.9 


0.3 


0.9 


0.3 


0.4_ 






_0.9_ 



Step 1: Matrices A and B are as follows: 

- 0.4 2 0.8 0.8 1 2 - 

1.25 2.5 2.5 1.25 1 1 

3 0.375 1 0.5 1 1.5 

2 1.5 1.2 0.75 0.75 1 

1 1.5 0.75 0.6 0.75 1 

5. 



Step 1: 



-1.60 


1.14 


1.23 


0.93 


1.26" 


1.42 


1.37 


0.87 


1.22 


3.52 


1.34 


0.90 


1.28 


2.23 


1.05 


.0.92 


1.08 


2.01 


78.5 


0.96. 



Step 2: 

Step 3: 

x = (0.92 0.90 0.87 0.93 0.96) 



L 1.2 1 2 2 1.2 1 



1 


1.5 


0.75 


0.6 


0.6 


1 


1.2 


1 


2 


2 


1.2 


1.5 


6 


3 


2 


1.5 


1.2 


1 


1 


1.2 


1.5 


2 


3 


6 


1 


1 


1 


1.33 


1.33 


1.33 


1.4 


1.4 


1.16 


1.16 


1 


1 


1 


1 


0.8 


0.8 


0.8 


1 


3 


1 


3 


1 


3 


2.25 



Step 4: 
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Step 2: 

Step 3: 

x=(l 1 0.75 0.6 0.6 1) 

Step 4: 

Step 5: 



A = 



0.4 0 
0 0 



0 0 
0 1 



0.5 0 0 

0 0 1 

0 0.75 0.6 0 1 

1 0 0 0 0 



Step 6: 

Step 7: In according to this step x(e 0 ) 5 = 0 and A is converted to 
"0.4 0 0 0 0" 



A = 



0 

0 

0 

0.75 

0 



0 

0.5 

0 

0.6 

0 



We notice that 5 $ R + , hence we can go to next step. 

Step 8: 

Step 9: J\ ew = {l}Jf w = { 6 }J 3 ew = {4 },J n 4 ew = { 6 }J n 5 ew = {1,3, 

4,6} and J" ew = { 2}. By pairwise comparison between four 
vectors (1 1 0 0.5 0 1), (0.4 1 0.75 0.5 0 1), 

(0.4 1 0 0.6 0 1) and (0.4 1 0 0.5 0 1), we gain 

x(e 0 ) = (0.4 1 0 0.5 0 1) as an optimal solution of Prob- 

lem (4a). 

Step 10: The optimal solution of the problem is x* = (0.4 1 
0.750.50.61) and the minimum value of the objective function 

| ? _ (0.4)(1)(0.5)(1) _ Q AAA 

^ ^ ~ (0.75)(0.6) — U 



5. Conclusion 

In this paper, we studied the monomial geometric program- 
ming problem with fuzzy relational inequality constraints de- 
fined by the max-product operator. Since the difficulty of this 
problem is finding the minimal solutions optimizing the same 
problem with the objective function Yl jGR +xJ\ we presented an 
algorithm together with some simplification operations to 
accelerate the problem resolution. At last, we gave three 
numerical examples to illustrate the proposed algorithm. 
Example (2) shows superiority of the presented approach with 
respect to preceding methods. 
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